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Abstract 

We apply an integral inequality to obtain a rigorous a priori estimate of the ac- 
curacy of the partial sum to the power series solution of the Riccati-Bernoulli 
differential equation. 

1 Introduction 

The deservedly celebrated Riccati-Bernoulli differential equation: 



/ 2 9 

y = X +y 



has the general solution [T]: 
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where c is an arbitrary constant and where 



J„{x) :-- 



X" 



1 - 



X 



+ 



22.l!-(n + l) 24 ■ 2! ■ (n + l)(n + 2) 



is the Bessel function of the first kind of order n, where n is any real number, and r(n + 1) is 
the gamma function. The very interesting history of (11. ip is detailed in Watson's standard 
treatise [5]. 

The equation f ll.ip is an example of a simple differential equation whose solutions form 
a family of transcendental functions which are essentially distinct from the elementary tran- 
scendents. 
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Unfortunately, the general solution (11. 2p does not easily lend itself to a rigorous error 
analysis of its accuracy in a particular interval of the independent variable. 

We will show how a simple application of an integral inequality allows one to estimate 
the accuracy of the partial sum of the Taylor series expansion of the solution within the 
latter's interval of convergence. 



2 Cauchy's theorem 

(In this section we follow [2j, Section IV.5). The general Cauchy Problem is to solve the 
ordinary differential equation initial-value problem: 



y' = f{x,y), y{xo):=yQ 



Let f{x,y) be expanded in the series 

/(a;, y) = ^ - xoYiy - yoY 



convergent for 



\x-xo\<Ri, \y-yo\<R2, {Ri > 0, i?2 > 0). 



(2.1) 



(2.2) 



(2.3) 



Then, according to Cauchy's theorem in the theory of differential equations, the problem 
(12. ip has a solution y{x) represented by the series 



(2.4) 



convergent in some neighborhood of the point xq. 

Cauchy's theorem allows one also to indicate the neighborhood of the point xq, in which 
the series (12. 4p converges. Namely, let M be a constant such that 




\f{x,y)\^M 



with 



|a; - a;o| ^ ri < \y - yo\ ^ r2 < R2, 

where i?i and R2 are numbers, defining the region (12.31) of the convergence of the series (12. 2p . 
and ri and r2 are some positive numbers. Then the series i \2A\) converges for 



\x — Xo\ < r. 



where 



r := ri(l 



-r2/2Mrl^ 



(2.5) 



It should be noted that the true interval of convergence is usually much larger than (12. 5p . 
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3 An integral inequality 

We use the notation of (12. ip . 

Theorem 1. Let I be the interval Xq ^ x ^ Xi, and suppose that f{x,y) > for all x E I , 
and that the differential inequality 



y' < f{xi,y{x)) 
also holds there. Then, the integral inequality 



f{xi,y{t)) dt 



'^4\dt^x-x, 



(3.i: 



(3.2) 



holds for all X E I . 

Proof. This is a simple consequence of the elementary calculus sufficient condition that a 
function be decreasing in an interval. Define 



g{x) :-- 



^}dt-{x-Xo). 



f{xi,y{t)) dt 

Taking the derivative and using the fundamental theorem of calculus, we obtain 

1 



9 [X) 



□ 



f{xi,y{x)) dx 

But, the inequality (13. ip shows that g'{x) ^ for all x E I. 
4 The Riccati— Bernoulli Initial Value Problem 

(This section is suggested by [2j, Section IV.5, but we add some important refinements.) 

Problem. It is required to find the first 11 terms of the power series expansion of the solution 
of the initial value problem 



I 2 , y 

y = X + — 



1/(0) := 



-1 



(4.1) 



and an interval of convergence for it. 



If we specialize the general solution (II. 2p to this initial value problem we find the following 
formula for the exact solution: 



y{x) 



X 



16 v^r(i)Ji(ix2) 



(4.2) 
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Unfortunately, the computation of the power series solution on the basis of the quotient 
fl4.2p . although theoretically possible, is computationally formidable. 

It is easier to use the equation f l4.ip to compute the derivatives of y{x) at x = directly: 

y' = x' + y^ = \, (4.3) 
y" = 2x + ^yy' = -^, (4.4) 

y'" = 2 + ly'' + \yy" = | , (4.5) 

y^'^ = \y'y" + \yy"' = -f^. (4.6) 

y^' = \y"' + yy"' + \yy^'^ = '^^^, (4.7) 

y^'^ = ^y"y"' + \y'y^'' + \yy''^ = -^^, (4.8) 



y''' = ^yV' + ^y"y^'' + \y'y^' + \yy''^ = , (4.io) 

= |(,(^))^ + 28y"'y^'^ + Uy"y('^ + Ay'y^'^ + \yy^'^ = ^ , (4.11) 

Therefore, by the formula (12. 4p . the first 11 terms of the series solution are: 

, 1 1 2 67 3 35 4 69 . 77 6 5317 . 

yix) = — 1 H — X X H a; x H x x H x 

^ 4 16 192 768 5120 15360 1032192 

2655 o 80641 „ 77851 ^ ^ 

x J x X + ■ • • f4 13) 

1835008 165150720 62914560 ^ ' ^ 

To find an interval of convergence of this series (14.130 . we use the Cauchy theorem. If 
X and y satisfy the inequalities 



then we may conclude that 



|x| ^ 0.5, + IK 1> 



\f{x,y)\ ^ |0.25[(y + l)-l]2 + x2| 
^ 0.25(|y + 1| + If + 
^ 1.25. 
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Therefore, in the formula (12. 5p . we may take 

ri:=0.5, r2 := 1, M := 1.25, 
and the value of r we obtain is: 

r = 0.5(1 - e"°-^) = 0.2753355 . . . 
Therefore, the power series solution ( I4.13P most certainly converges for \x\ ^ 0.27. 

5 The Accuracy of a Partial Sum from the Integral Inequality 

We will consider the following concrete problem, although the principles are of general ap- 
plicability. 

Problem. It is required to determine the accuracy of the partial sum of degree 9 of the power 
series solution (I4.13P in the interval ^ x ^ 0.2. 

Since the series (14.131) is the Maclaurin expansion of y{x), we must estimate the re- 
mainder term RqIx), which we write in the Lagrange form: 



where < Gg < 0.2. We have to estimate ?/^^°^(x), the formula for which is given in (I4.12p . 
for all values of x in the interval ^ x ^ 0.2. The formulas (14. 3 p through (14. lip show us, 
finally, that we must estimate y{x) itself in ^ x ^ 0.2. 

The estimate of y{x) via our integral inequality (13. 2 p constitutes the novelty in this paper. 

Maintaining the notation of (13. 2p . we see that 

Xq := 0, Xi := 0.2, 

that the right hand side 

f{x,y):=x' + !^>0 
on /, and the differential inequality (13. ip becomes: 

^^0.04+^ 
dx 4 

Therefore, the integral inequality (13. 2 p becomes 

^dt<x. 



l0.04+Wt)2 
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But 



Therefore, 



or 



dy/dt 



0.04 + \y{tY dt 



d 

— < 10 arctan 



0.4 



yix) (-1) 
10 arctan 10 arctan ^ x. 



arctan 



0.4 

yit) 



0.4 



X (-1) 

— - ^ h arctan — , 

0.4 10 0.4 ' 



and taking the tangent of both sides and reducing, we obtain the estimate 




which holds for all x E I. 

The function on the right-hand side of (15.11) is monotonically increasing in /, and 



tan(f ) - 1 



1 + 



0.2 



-0.9447608 ■ ■ ■ < -0.94, 



and we have proved that the following inequality is true for all x E I: 



-1 < y{x) ^ -0.94 



{Note: the true value of y{0.2) is 



2/(0.2) 



-0.9497771 



so the estimate (15. ip . with an error of —0.00501 . . . , or about 0.53%, is quite good!) 



Now we must estimate y'{x) 



(10) I 



x) using the formulas (14. 3 P through (I4.12p . 



Theorem 2. The following inequalities are valid for all x E I: 



-1 




y{x) 




-0.94, 


0.22 




y\x) 




0.29, 


-0.15 




y"{x) 




0.3, 


1.87 




y"\x) 




2.12 


-1.13 




l/(^)(x) 




-0.74, 


1.17 




y^^\x) 




1.93, 


-3.38 




y(6)(a;) 




2.23, 


14.59 




y(^)(x) 




27.12, 


-61.96 




y(«)(a;) 




-22.73 


92.03 




y(9)(a;) 




146.76, 
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and finally 



-665.9 < y^^^\x) < 281 



Proof. All of our estimates come from worst-case values applied to each of the summands in 
the formulas. 
By gJD, 

y =x + — . 

Therefore, using (15. 2p . we conclude that for all x G /, 



y' < 0.2^ 



4 



0.29, 



while 



> 0^ + ^ ' ' = 0.2209 > 0.22. 



Therefore, we obtain the bounds 



0.22 < y'{x) < 0.29 



By gaD, 

,, 1 , 

y =2x + -yy . 

Therefore, using (15. 3p . we conclude that for all x G /, 

y" < 2(0.02) + -(-0.94)(0.22) = 0.2966 < 0.3, 
2 

while 



y" > 2(0) + -(-1)(0.29) = -0.145 > -0.15. 



Therefore, we obtain the bounds 



(5.3) 



-0.15 < y"{x) < 0.3 



By 



III , ^ 12 , ^ II 

y =2+2^ ^ ■ 



Therefore, using (15. 3p and (15. 4p . we conclude that for all x G /, 

y'" < 2 + ^(0.29)2 + ^(-0.15)(-1) = 2.11705 < 2.12, 



(5.4) 
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while 



y'" > 2 + ^(0.22)2 + i(0.3)(-l) = 1.8742 > 1.87. 



Therefore, we obtain the bounds 



By dUD, 



1.87 < y"{x) < 2.12 



- ^yy + ^yy 



Therefore, using (15. 3p . (15. 4p and (15. 5p . we conclude that for all x E I, 

y^^^ < ^(0.29)(0.3) + i(-0.94)(1.87) = -0.7484 < -.074, 



while 



> ^(0.29)(-0.15) + i(-l)(2.12) = -1.12525 > -1.13. 



Therefore, we obtain the bounds 



-1.13 < y^^\x) < -0.74 



By (HZ!), 



y^ - ^y +^yy +^yy^' 



Therefore, using (15.31) through (15.61) . we conclude that for all x E I, 

y^^^ < ^(0.3)2 + 2(0.29)(2.12) + ^(-1)(-1.13) = 1.9296 < 1.93, 



while 



3 1 

^(5) > _(o)2 + 2(0.22)(1.87) + -(-0.94)(-0.74) = 1.1706 > 1.17. 
2 2 



Therefore, we obtain the bounds 



1.17 < y^^\x) < 1.93 



By dH), 



= 5y"y"' + \y'y^'^ + \yy^''^ 



Therefore, using (15. 3p through (15. 7p . we conclude that for all x 
y'^^^ < 5(0.3)(2.12) + ^(0.22)(-0.74) + ^(-0.94)(1.17) = 

while 

> 5(-0.15)(2.12) + ^(0.29)(-1.13) + ^(-1)(1.93) = 
Therefore, we obtain the bounds 

-3.38 < < 2.23 . 

By (USD, 

Therefore, using (15. 3p through (15. Sp . we conclude that for all x E I, 

y^^^ < 5(2.12)2 ^ ^(_o.l5)(-1.13) + 3(0.29)(1.93) + ^(-l)(-3.38) 
= 27.11235 < 27.12, 

while 

y(^) > 5(1.87)2 + 1^(0.3)(-1.13) + 3(0.22)(1.17) + ^(-1)(2.23) 
= 14.5992 > 14.59. 

Therefore, we obtain the bounds 

14.59 < y^'^\x) < 27.12 . 

By dHO]), 

Therefore, using (15. 3p through (15. 9p . we conclude that for all x E I, 

35 21 7 1 

^(8) < _(i.87)(-0.74) + y (0.3)(1.93) + -(0.29)(2.23) + -(-0.94)(14.59), 

= -22.73085 < -22.73, 

while 

> ^(2.12)(-1.13) + y (-0.15)(1.93) + ^(0.29)(-3.38) + ^(-1)(27.12) 
= -61.95345 > -61.96. 



= 2.2231 < 2.23, 

= -3.37425 > -3.38. 
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Therefore, we obtain the bounds 



61.96 < y^^\x) < -22.73 



(5.10) 



By gUD, 



y 



(9) 



35 



Therefore, using (15. 3p through (15.101) . we conclude that for all x G / 



yy 



(8) 



^(9) ^ rr (_i.i3)^ + 28(2.12)(1.93) + 14(0.3)(2.23) + 4(0.29)(27.12) + ^(-1)(-61.96) 



35 



while 



^(9) ^ rr (_o.74)^ + 28(1.87)(1.17) + 14(-0.15)(2.23) + 4(.22)(14.59) + ^(-0.94)(-27.73) 



146.75575 < 146.76, 



35 



= 92.0335 > 92.03. 
Therefore, we obtain the bounds 



92.03 < y^^\x) < 146.76 



(5.11) 



By gl2D, 

^(10) = 63^(4)^(5) + ^^y.yie) ^ 18^.^(7) + 9^,^(8) + 1 ^^(9) 
Therefore, using (15.31) through (15.111) . we conclude that for all x G /, 



< 63(-0.74)(1.17) + 42(2.12)(2.23) + 18(0.3)(27.12) + ^(0.22)(-22.73) 

+ ^(-0.94)(92.03) 
= 280.9922 < 281, 



while 



y^^^^ > 63(-1.13)(1.93) + 42(2.12)(-3.38) + 18(-0.15)(27.12) + ^(0.29)(-61.96) 

+ ^(-l)(146.76) 
= -665.8137 > -665.9. 
Therefore, we obtain the bounds 
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-665.9 < y^^^\x) < 281 



The inequahty (I5.12p was the goal of this long and detailed computation! 



(5.12) 



□ 
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Now we can state the accuracy of the partial sum. 



Theorem 3. For all x in ^ x ^ 0.2, the partial sum of degree 9: 



2655 o 

x^ + 



67 o 

x'^ 

192 

80641 



35 
768 



69 
5120 



77 



■ X 



15360 



5317 
1032192 



X 



1835008 



165150720 



X 



(5.13) 



approximates the true value, y{x), of the solution series (14.131) . with an error that does not 
exceed 2 units in the eleventh decimal place. 



Proof. The estimate 
\Rg{x)\ = 
completes the proof. 



10! 



X 



10 



^^(0.2)1° 
10! ^ ' 



1.878...- 10-11 < 2 - 10 



-11 



□ 



6 Conclusions 

Our integral inequality (13. 2p can be applied to wide classes of differential equations. 

For example, our method allows us to prove that in the interval ^ x ^ 0.4 the polyno- 
mial 

/\ ^9 '^Q 1/1 

y(x) := IH \ X -\ x^ -\ x -\ x^ 

' 4 16 192 96 200 

approximates the true solution, y{x), of the initial value problem 

4y' = x + y\ 2/(0) :=1, 

with an error that does not exceed 2 units in the fifth decimal place. 

We chose the Riccati-Bernoulli equation because it illustrates the process so perfectly 
and because a direct estimate of the accuracy of the partial sum of the series solution is 
troublesome. 

We did not investigate the accuracy for the negative half of the interval, i.e., for —0.2 ^ 
x ^ 0, which we leave as an exercise for the reader. The only change occurs in the estimate 
of y"{x) since then x can take negative values. 

Finally, we observe that we did not exploit the sign of the error. In fact, our estimate 
(I5.12P allows us to say that the error we commit is between an error in defect smaller than 
2 units in the eleventh decimal place and an error in excess smaller than 8 units in the twelfth 
decimal place. Therefore we can centralize the error by adding the term 

281 - 665.9 10 _ 1283 
2 ■ 10! ^ ~ "24192000^ 

to our polynomial (I5.13P and obtain an approximating polynomial whose maximum error is 
between ±1.4 units in the 11th decimal place. 
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